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MHD-MODELS OF QUIESCENT PROMINENCES 
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^ • I ntroduction . it is well known the close connection 
between the sudden disappearance ( "disparition brusque”) of 
the quiescent prominences and the two— ribbon flares (see 
e-g. [6]). During this dynamic phase the prominence ascends 
rapidly (typically with a velocity about 100 Km/sec) and 
disappears. in another later stage is observed material 
falling back into the chromosphere. The impact of this 
downfalling matter on the chromosphere produces the 
Ho<, -brigthening , which shows the symmetric double pattern. 

The occurrence of the "disparition brusque" is thought to be 
a consequence of a plasma instability of magnetohydrostatic 
structures (see e.g. [11]). Two-dimensional MHD-models for 
quiescent prominences have been worked out since the fifties. 
They describe the prominence in magnetohydrostatic 
equilibrium under the action of Lorentz forces, gas pressure 
gradients and the gravitational force. However, the 
stability properties of most of these models are not yet 
determined. We analyze by means of the MHD-energy principle 
[2] the stability properties of four prominence models. We 
show that all considered models undergo instabilities for 
parameters outside of the observed range at quiescent 
prominences. We consider the possibility that such 
instabilities in the flare parameter range may indicate iust 
the onset of a flare. 

— — ^$ u ^-j- Dr -*- uiri an( ^ S t ability Theory, vie define a coordinate 
system with x-axis along prominence, y-axis perpendicular to 
the prominence sheet and z-axis vertical (opposed direction 
of the gravity acceleration) . We take into account only 
plasma structures, whibh are independent on x. In a 
two-dimensional theory the magnetic field can be expressed 
as : 

B=VACy.Z)xg x + B >t (y,z) e„ (1) 

where A is the x-component of a vector potential. The equili- 
brium condition, which A must satisfy reads: 

AA = -4ir|^M) , n(A»=PCA,^)4-i f B*(A) (2) 

where P is the plasma pressure, the external gravitational 
field and J =- ‘On/S^s-'DP/d^ the mass density. Any 
two-dimensional prominence model corresponds to a particular 
choice of the functions P(A,0), Bx (A) and the boundary 
conditions. In order to analyze the stability properties of 
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prominence models we use the MHD-energy principle of 
Bernstein et al . (1958), according to which the stability 

of an equilibrium configuration is determined by the 
behaviour of the potential energy functional ) 

resulting from a perturbation §(r,t) to the system. For the 
two-dimensional equilibrium class the functional 6 "W {%,£*) 
may be written m the form: 

+ dev,. Im^ 

&] } dr 131 

where 7 x= e y f- + & V X A , E x = VxA*g* . We have 

supposed complex three-dimensional displacements of the 
following generalized form: 

iwt r- _ lkx + iwt 

g(r,t)=£(r)6 =Q^ x (y, z)e x 4-g A (y,z)] e (4) 

It is assumed periodic boundary conditions along the promi- 
nence axis (x-direction) and % =0 on the edge of the plasma 
region in the y,z-plane. The x-integration in equation (3) 
is to be carry out over one period. If 8 W is positive for 
all displacements v/hich satisfy the boundary conditions, 
then the equilibrium is stable. The MHD-energy principle is 
a necessary and sufficient criterion for stability. The 
stability problem reduces therefore to analize the sign of 
the minimum of Sw. It is interesting to note that to the 
variational problem which implies the minimisation of 5w, 
is associated the Euler-Lagrange equation: 

- 5 co | , cr)=P(?(r)) (5) 

where F is a self-adjoint differential operator v/ith time 
independent coefficients (see e.g. [1]). We use the 
normalization constraint: h 5 $ lf|| 2 d 3 r * 1 and obtain the 

spectrum of eigenvalues further one has co^, 1(1 = Mm<$W(g,S*) 


3. Stability Results for Prominence Models. We have 
evaluated the energy principle for four prominence models: 
Menzel, (M) , [10]; Dungey, (D) , [3]; Kippenhahn and 

Schliiter , (KS) , [7]; Lerche and Low, (LL) , [9]. One obtains 

generally stability statements by two procedures: (a) 

analytically, by manipulating the energy functional <$W to 
recognize a definite sign and so to infer about the 
stability properties, (b) numerically, by carring out the 
minimisation of <Sw with the aid of a computer code. Our 
analytical results concerning to these models are reported 
elsewhere [4] . We obtain global stability for the KS-model 
in case of arbitrary 3D-displacements [5] ; for the other 
models the stability statements are restrictive to special 
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classes of displacements (2D, long wavelengths, etc.). In 

order to obtain more general results concerning to broader 

class of displacements we have developed a numerical code 

based on the finite element-method. This procedure provides, 

besides qualitative stability statements, the largest 

growth rate in case of instability and the frequency of the 

fundamental oscillations of stable systems. Our code was 

tested successfully by applying it to simple systems, whose 

dynamic properties are well known (e.g. AlfvSn waves, 

current sheets and sound waves in a constant gravitational 

field) . We obtain numerically a stable behaviour for the 

four models in the parameter range of the observations of 

typical quiescent prominences (T=7.10 3o K, n =5 . 10 1 0 cm“ 3 , 

r e r 

B=5G, thickeness Y=5.10 3 Km, height Z=1.5-5.10 4 Km). All 
models describe horizontal large-scale oscillations with 
periods between 16 and 80 min (see Table 1). Reported data 
indicate that quiescent prominences undergo actually 
horizontal oscillations from and towards the perturbation, 
which is originated in solar flares. The observed periods 
range from 6 to 80 min. [2], [8], [1]. 


Table 1. Periods of horizontal 
oscillations in quiescent 
prominences . 


Model 


LL model 


Period (min) 
40 

55-80 

16 

17-50 


Unstable behaviour is found 
only out of the observed 
parameter range. When one 
considers typical 
parameters of flares (e.g. 
T~10 7o K, n ~ lO^cm" 3 , 
B~50 G) oni obtains 
A s maXL > A . Just for 

f -n2 ai 


17-50 these parameters we get 

instability only in case of 
The other models continue stable around these 


parameters. Figure 1 shows function of eks/i/Z (we use 
dimensionless variables, so that CO* . is normalized by g/h, 

where g=2 . 74x10 4 cm/sec 2 , h=kT/mg is the density scale height 
and m is the proton mass) • We have separated the different 
physical effects (electromagnetic, compressional and 
gravitational parts) in the energy functional 

so that it is possible to infer about the 


nature of the instability. For the situation illustrated on 
Fig. 1 we find that the instability is driven mainly 
by electromagnetic forces. Gravitation provides also an 
instabilizing effect. In opposition to this, the 
compression has continually a stabilizing effect. Such a 
gravitational-electroma gnet ic mixed instability has a 
typical growth rate r = Jf/h 1 |w mw | /2lf ~ 4.77 \Q~\ec l , i- e. 4 
growth timet'** 5.8 h. On the other hand, the impulsive phase 
of a flare elapses within few minutes, so that the found 
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growth time may explain rather the 
evolution of the flare ribbons and the 
associated loop system which last 
several hours . The cause of the 
parameter shift in the quiescent 
prominence may reside in an external 
perturbation in form of a shock wave 
generated by a distant flare or in an 
internal perturbation in form of a 
newly emerging flux in the same active 
region. Our results can be considered 
as preliminaries because the studied 
two-dimensional models are still very 
simple to describe the complexity 
of quiescent prominences. However, 
further applications of our stability 
method to more realistic models' give the possibility of a 
better description of the impulsive phase of solar tiares. 



Fig . 1 . LL-Model : 

minimum eigenvalue 

G J 2 . as function 
min 

of <*=/3/ 2. 
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